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ABSTRACT 

We present a search for O I in the spectra of nine 4.9 < zqso < 6.4 QSOs 
taken with Keck/HIRES. We detect six systems with Nq, i > lO^'^'^ cm~^ in the 
redshift intervals where O I A1302 falls redward of the Lya forest. Four of these 
lie towards SDSS J1148+5251 (;zqso = 6.42). This imbalance is unlikely to arise 
from variations in sensitivity among our data or from a statistical fluctuation. 
The excess O I occurs over a redshift interval that also contains transmission in 
Lya and Ly/5. Therefore, if these O I systems represent pockets of neutral gas, 
then they must occur within or near regions of the IGM that are highly ionized. 
In contrast, no O I is detected towards SDSS J1030+0524 (zqso = 6.30), whose 
spectrum shows complete absorption in Lya and Ly/5 over I\z ~ 0.2. Assuming 
no ionization corrections, we measure mean abundance ratios ([0/Si]) = — 0.04± 
0.06, ([C/0]) = -0.31 ± 0.09, and ([C/Si]) = -0.34 ± 0.07 (2(7), which are 
consistent with enrichment dominated by Type II supernovae. The 0/Si ratio 
limits the fraction of silicon in these systems contributed by metal-free very 
massive stars to < 30%, a result which is insensitive to ionization corrections. 
The ionic comoving mass densities along the -Zqso > 6.2 sightlines, including only 
the detected systems, are = (7.0 ± 0.6) x 10-^ Qsm = (9.6 ± 0.9) x 10"^ 
and fie II = (1.5 ±0.2) x 10"^ 
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1. Introduction 

The state of the intergalactic medium (IGM) at redshift 2; ~ 6 remains under consider- 
able debate. Significant transmitted flux in the Lja forest at 2; ~ 5.7 means that the IGM 
must have been highly ionized by at least ~ 1 Gyr after the Big Bang (Becker et al. 2001; 
Djorgovski et al. 2001; Fan et al. 2001, 2004). Each of the four known QSOs at ^qso > 6.2 
show Gunn-Peterson troughs (Gunn & Peterson 1965) over at least a narrow redshift interval 
(Becker et al. 2001; Pentericci et al. 2002; Fan et al. 2001, 2003, 2004; White et al. 2003). 
This complete lack of transmitted flux has been interpreted as an indication that the tail end 
of cosmic reionization may extend to z ~ 6. However, Songaila (2004) found the evolution 
of transmitted flux over 2 < z < 6.3 to be consistent with a smoothly decreasing ionization 
rate and not indicative of a sudden jump in the Lya optical depth at z ^ 6. 

Significant variations in the fraction of transmitted flux are common among sightlines at 
the same redshift (Songaila 2004). While the spectrum of SDSS J1030+0524 (^qso = 6.30) 
shows complete absorption in Lya and Ly/? over a redshift interval Az ^ 0.2 (White et al. 
2003), transmitted flux appears over the same redshifts in the Lya, Ljf3, and Ly7 forests 
of SDSS J1148+5251 {zqso = 6.42) (White et al. 2005; Oh & Furlanetto 2005). Either 
the IGM is highly ionized everywhere at 2; ~ 6.3 and long stretches of complete absorption 
are the result of fine blending, or the neutral fraction of the IGM is patchy on large scales. 
A patchy IGM could result from the dumpiness of the IGM and the clustering of ionizing 
sources (Furlanetto & Oh 2005). 

Studies of H I transmitted fiux are ultimately hampered by the large optical depths 
expected for even a small neutral fraction. The presence of a Gunn-Peterson trough can at 
best constrain the volume- and mass- weighted H I neutral fractions to > 10~^ and > 10~^, 
respectively (Fan et al. 2002). Alternative measurements are required to probe larger neutral 
fractions. The non-evolution in the luminosity function of Lya-emitting galaxies provides an 
independent indication that the IGM at 2; ~ 6.5 is highly ionized (e.g., Malhotra & Rhoads 
2004; Stern et al. 2005), although Lya photons may escape as a result of galactic winds 
(Santos 2004) or from locally ionized bubbles created by clustered sources (Furlanetto et al. 
2004a). In the future, observations of redshifted 21 cm emission/absorption should trace the 
growth and evolution of ionized regions at 2; > 6, placing strong constraints on reionization 
scenarios (e.g., Tozzi et al. 2000; Kassim et al. 2004; CariUi et al. 2004; Furlanetto et al. 
2004b). 

A high neutral fraction must also have a measurable effect on metal absorption lines (Oh 
2002; Furlanetto & Loeb 2003). Overdense regions of the IGM should be the first to become 
enriched, due to the presence of star-forming sources, yet the last to remain ionized, due to the 
short recombination times (Oh 2002). Low- ionization metal species should therefore produce 
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numerous absorption features in the spectra of background objects prior to reionization. A 
particularly good candidate is O I, which has an ionization potential nearly identical to that 
of H I and a transition at 1302 A that can be observed rcdward of the Lya forest. Oxygen 
and hydrogen will lock in charge-exchange equilibrium (Osterbrock 1989), which ensures 
that their neutral fractions will remain nearly equal, 

^Hi noi /.v 

/hi = ~ . (1 

nn no 

Despite the increased photo-ionization cross section of O I at higher energies, this relationship 
is expected to hold over a wide range in /hi (Oh et al. 2005, in prep). 

In this work, we present a search for O I in the spectra of nine 4.9 < zqso < 6.4 
QSOs. This is the first time a set of zqso > 5 spectra has been taken at high rcsohition 
(R = 45,000). Our sample includes three objects at -Zqso > 6.2, where wc might expect to 
see an "O I forest" (Oh 2002). In §2 we describe the observations and data reduction. The 
results of the O I search are detailed in §3. In §4 we demonstrate a significant overabundance 
of O I systems towards the highest-redshift object, SDSS J1148-I-5251 {zqso — 6.42), and 
compare the overabundance to the number density of lower-redshift O I systems and other 
absorbers with high H I column densities. Measurements of the relative metal abundances 
for all the detected systems are described in §5. In §6 we discuss the significance of these 
excess systems for the enrichment and ionization state of the IGM at z ~ 6. Our results 
are summarized in §7. Throughout this paper we assume fl = 1, fl^i = 0.3, f^A = 0.7, and 
Ho = 70 km s"^ Mpc'^ 



2. The Data 

Observations using the Keck HIRES spectrograph (Vogt et al. 1994) were made between 
2003 February and 2005 June, with the bulk of the data acquired during 2005 January and 
February. Our QSO sample and the observations are summarized in Table 1. All except the 
2003 February data were taken with the upgraded HIRES detector. We used an 0'.'86 slit, 
which gives a velocity resolution FWHM of Av = 6.7 km s~^. 

The continuum luminosities in our sample approach the practical detection limit of 
HIRES. In addition, the spectral regions of interest for this work lie in the far red, where 
skyline contamination and atmospheric absorption are major concerns. To address these 
difficulties, the data were reduced using a custom set of IDL routines written by one of 
us (GDB). This package uses optimal sky subtraction (Kelson 2003) to achieve Poisson- 
limited residuals in the two-dimensional sky-subtracted frame. One-dimensional spectra 
were extracted using optimal extraction (Horne 1986). Blueward of the Lya emission line. 
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the continuum for these objects is highly obscured due to strong absorption from the Lya 

forest. Furthermore, single exposures were typically too faint to fit reliable continua even in 
unabsorbed regions. Individual orders were therefore combined by first fiux calibrating each 
order using the response function derived from a standard star. The combined spectra were 
continuum-fit by hand using a cubic spline. Standards stars were also to correct for telluric 
absorption. Estimated residuals from continuum fitting (redward of Lya) and atmospheric 
correction are typically within the flux uncertainties. 

3. O I Search 
3.1. Technique 

At 2; > 5 it becomes increasingly difficult to identify individual Lya lines due to the 
large opacity in the Lyct forest. Therefore, candidate O I lines must be confirmed using other 
ions at the same redshift. Si ll A1260 and C ll A1334 are the most useful for this purpose 
since their rest wavelengths place them redward of the Lya forest, at least over a limited 
range in redshift. These transitions are also expected to be relatively strong over a range 
in ionization conditions. At lower redshifts where Si II A1260 has entered the forest, we are 
still sensitive to Si ll A1304, although it is weaker than Si ll A1260 by a factor of 7. 

We searched for O I in each sightline over the redshift range between the quasar redshift 
Zqso and the redshift z„^i^ = (1 + 2;qso)Aq;/Ao i — 1, where O I enters the Lya forest. To 
identify a system, we required absorption at the same redshift in O I and either Si II A 1260 or 
C II A1334, and that the line profiles for at least two of the ions be comparable. In practice. 
Si II A1260 proved to be the best indicator over the lower-redshift end of each sightline, since 
the signal-to-noise {S/N) ratio was typically better over that region of the spectrum. At 
redshifts where Si ll A1260 was lost in the forest we relied solely upon C ll for confirmation. 
By chance, the two O I systems we identified using only C ll were also strong enough to be 
detected in Si 11 A1304. 

We identified six O I systems in total along our nine fines of sight. Their properties 
are summarized in Table 2. Four of the six systems appear towards SDSS J1148-I-5251, and 
we examine this potential overabundance in detail in §4. Voigt profiles were fitted using 
the VPFIT package written, by R. Carswell. Due to line blending and the modest S/N 
ratio over some regions of the data, it was often necessary to tie the redshift and/or Doppler 
parameter of one ion to those of another whose profile was more clearly defined in order to 
obtain a satisfactory fit. In practice, the resulting total column densities are very similar 
to those obtained when all parameters are allowed to vary freely, as would be expected for 
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systems still on the linear part of the curve of growth. Comments on individual systems are 
given below. 

3.1.1. SDSS J0231-0728: ^sys = 5.3364 (Figure 1) 

This system shows a central complex in O I, Si ll, and C ll that is best fit using three 
components with a total velocity span of Ai> 38 km s~^. We measure the redshifts and h- 
parameters of the bluest two central components from Si ll. The redshift and 6-parameter for 

the reddest central component are measured from O I, which is least affected by additional 
lines towards the red. Leaving all parameters untied produces a two-component fit for 
C II, which is likely due to low S/N . However, the total column densities remain virtually 
unchanged. 

A single component at Aw ^ —115 km s~^ appears in Si ll and C ll but not in O I. 
Since O I is more easily ionized than either Si ll or C ll {AEq i = 13.6 eV, Ai^si n = 16.3 eV, 
AEq 11 = 24.4 eV) this component probably arises from gas that is more highly ionized than 
the gas producing the O I absorption. Three additional components appear redward of the 
central O I complex. These are visible in Si ll, C ll, and Si IV and presumably correspond to 
gas in a still higher ionization state. We use Si ll to constrain the redshift and 6-parameters 
for all the non-0 I components, although fitting without tying any parameters produces very 
similar results. 

3.1.2. SDSS J1623+3112: z^ys = 5.8408 (Figure 2) 

This system shows multiple components with a total velocity span of Av a; 81 km s~^. 
The Si II is strong enough that the 1304 A transition is clearly visible. Some additional 
absorption occurs along the blue edge of the O I complex. However, since no corresponding 
absorption appears in C ll this features is unhkely to be genuine O I. We use C ll to fix 
the redshifts and 6-parameters of the bluest two components and the reddest component for 
all ions. The z and b values for the remaining components are tied to O I. The weakest 
C II component was automatically dropped by the fitting program. However, this should 
not have a significant impact on the summed C ll column density. No Si IV was detected 
down to iVsiiv > 10^^ cm"^. 

A strong line appears at a velocity separation Av ~ —180 km s~^ from the O I complex. 
Weak absorption appears at the same velocity in Si II A1304. However, the corresponding 
C II region is completely covered by Mg I A2853 from a strong Mg ll system at z^hs = 2.1979. 
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Therefore, while these features are plausibly associated with the O I complex, their identity 
cannot presently be confirmed. We do not include these lines when computing the relative 
abundances, although their O l/Si ll ratio is consistent with those measured for the other 
O I systems. 



3.1.3. SDSS J1148+5251: z,y^ = 6.0097 (Ftgure 3) 

This system shows multiple components with a velocity span of Ai> 41 km s~^. The 
Si II A1260 transition falls in the Lya forest. However, Si ll A1304 is clearly visible. The 
C II occurs among moderate-strength telluric absorption. In correcting for the atmosphere 
we recovered a C II velocity profile very similar to that of O I and Si II. Some additional 
absorption around the C II arc probably telluric lines that were not completely removed. Wc 
achieved the best fit for O I using four components while allowing the continuum level to 
vary slightly. Using fewer components resulted in virtually the same total A^ou which is to 
be expected if the hues he along the hnear part of the curve of growth. We used the z and h 
values for these four O I lines to determine the corresponding column densities for Si ll and 
C II. Leaving all parameters untied resulted in small changes in the 6-parameters and in the 
preferred number of lines used by VPFIT (3 for Si II instead of 4), but negligible changes in 
the total column density for each ion. We detect no Si IV down to A^"si iv > 10^^'^ cm~^. This 
system docs show C IV absorption in the low-resolution NIRSPEC spectrum taken by Barth 
et al. (2003). However, the kinematic structure cannot be determined from their data, and 
it is therefore unknown whether the C IV absorption arises from the same gas that produces 
the lower-ionization lines. 



3.1.4. SDSS J1148+5251: Zsys = 6.1293 (Figure 4) 

The Si II A 1260 line for this system falls in the Lja forest. However, Si II A 1304 is 
strong enough to be visible. We tie the redshifts and 6-parameters for Si ll and C ll at 
z — 6.129267 to those for O I, which has the cleanest profile. Leaving all parameters untied 
produces changes in the 6-parameters that are within the errors and negligible changes in 
the column densities. Possible additional C ll lines appear at Av — —67, -29, and 26 km 
from the strong component. These do not appear in O I and would probably be too weak 
to appear in Si II A1304. Wc do not include them when computing the relative abundances 
for this system. No Si IV was detected down to A^siiv > 10^^'^ cm~^. 
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3.1.5. SDSS Jl 148+5251: Zsys = 6.1968 (Figure 5) 



The Si II A1260 and C II lines for this system occur in transparent gaps within regions 
of mild telluric absorption. In some of our exposures, the O I line falls on the blue wing 
of a strong telluric feature. However, the O i can be clearly seen in the uncorrected 2005 
January data, where the heliocentric velocity shifted the observed QSO spectrum further to 
the blue. We tie the redshifts and 6-values for O I and C ll to those of Si ll, for which the 
line profile is the cleanest. Leaving all parameters untied produces negligible changes in the 
column densities. 



We use a single component for the primary fit to this system. The C ll falls on a sky 
emission line and hence the line profile is rather noisy. We therefore tied the redshift and 
6-parameter for C ll to those of Si ll, while z and fo for O I were allowed to vary freely. The 
fitted z and h values for O I agree with those for Si ll within the fit errors. Tying z and h for 
O I to those of Si ll produced no significant change in A'^o i- Both Si ll A1260 and O I are well 
fit by a single line, while the blue edge of the C II shows a possible additional component. 
We fit this component separately but do include it when computing relative abundances. 

4. Overabundance of O i systems towards SDSS J1148+5251 

Despite the relatively small total number of detected O l systems, a noticable overabun- 
dance appears towards SDSS J1148+5251. It should be noted that this object has one of 
the brightest continuum magnitudes among QSOs at 2; > 5 (Fan et al. 2003). We also spent 
considerable integration time on this object. Therefore, the data are among the best in our 
sample. In this section we address the null hypothesis that the apparent imbalance in the 
distribution of O I systems results from a combination of incompleteness and small number 
statistics. 

For each sightline we used an automated scheme to determine the absorption pathlength 
interval AX over which we are sensitive to systems above a given O I column density. For a 
non-evolving population of sources expanding with the Hubble fiow, we expect to intersect 
the same number of sources per unit X at all redshifts, where X is defined for = 1 as 



3.1.6. SDSS J 1148+5251: z^y^ = 6.2555 (Figure 6) 



dX 




dz . 
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At the redshift of each pixel where we are potentially sensitive to O I, we inserted artificial 
single-component systems containing O I, Si ll, and C ll with relative abundances set to 
the mean measured values (see §5) and a turbulent Doppler width of 6 = 10 km s~^. The 
minimum Nq j at which the system could be detected by the automated program was then 
determined. A detection required that there be significant (4(t) absorption at the expected 
wavelengths of O I and at least one of the other two ions. The redshifts of at least two of 
the three fitted profiles must also agree within a tolerance of A.Ztoi — 0.0002, and the fine 
widths must agree to within a factor of two. In the case where the fitted profiles agree for 
only two of the ions, the spectrum at the wavelength of the third ion must be consistent with 
the expected level of absorption. This allows for one of the ions to be blended with another 
line or lost in Lja forest. Lines that fell in the Fraunhofer A band were not considered. 
This method detected all of the known O I systems and produced only one false detection 
in the real data, which was easily dismissed by visual inspection. We note that this method 
focuses on our ability to detect systems similar to those observed towards SDSS Jl 148+5251. 
Multiple-component systems such as the two at lower redshift are in general easy to identify, 
even in noisy regions of the spectrum. 

The results of the sensitivity measurements are shown in Figure 9. The solid lines show 
the AX over which we are sensitive to systems with A'^o i > A'o ™- The O I column densities 
for the four systems seen towards SDSS J1148-I-5251 are marked with circles. The dashed 
fine shows the mean AX as a function of A'q™ for all nine sightlines. While it is clear that 
the path length over which we are sensitive to low-column density systems is considerably 
smaller for some sighthnes than for SDSS J1148-I-5251, there does exist significant coverage 
across the total sample even at the A^o i of the weakest system in our sample. 

Completeness effects will depend on the intrinsic distribution of line strengths, in that 
weaker systems will more easily be lost in \ow-S/N sightlines. It is useful to define an effective 
absorption pathlength interval AX^^ such that, for a sightline with a total pathlength AX*°* 
and sensitivity S{N), we would expect to detect the same number A/" of systems with A^ > 
Ainin over AX^^ if we had perfect sensitivity. For a column-density distribution f{N) — 
dW/dNdX, this gives 



/■oo poo 

AX^^ / f{N) dN = AA:*°* / f{N) S{N) dN , 



(3) 



where 5'(A^) is the fraction of the total pathlength over which we are sensitive to systems 
with column density A^. The sensitivity function can be determined directly from Figure 9. 
We assume a conventional power-law distribution in O I column densities /(A^oi) c< A'q". 
Hence, 



.Jn N~^S{N)dN 



aa:"** = aa:*°* ^^^^^ — — . (4) 

N-^dN ^ ' 
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Taking into account our measured sensitivity, the maximum likelihood estimator for a for 
all six O I systems is a = 1.7. A K-S test allows 1.3 < a < 2.1 at the 95% confidence level. 

In Figure 10 we plot /S.X^^ as a function of a taking A^min = 10^^-^ cm~^ as the cutoff 
value. The solid hne shows AX^^ for SDSS J1148+5251, while the dashed hne shows /^X^^ 
for the remaining eight sighthnes, divided by 5.7 to match SDSS J1148+5251 at a = 1.7. 

It is easy to sec that the ratio ^X^^^/ ^X^-^^^ does not change significantly with a (< 4% 
over 1.1 < a < 3.0). This is purely a coincidence, most likely related to the fact that the 
overall high quality of the SDSS J1148+5251 data is balanced by the placement of the O I 
systems in a difficult region of the spectrum. However, it means that the relative numbers 
of expected O I detections for SDSS J1148+5251 and the remaining eight sightlines should 
be nearly invariant over a wide range in plausible column density distributions. 

We can now evaluate the likelihood that the imbalance of O I detection results from 
small number sampling. For a Poisson distribution, the likelihood of detecting M systems 
when the expected number of detections is v is given by 

P{Af,iy)^-^. (5) 

If we assume that the expected number of detections towards SDSS J1148+5251 is u = 
6 (AXlf^g/AXlfi) ~ 0.90, then the probability of detecting at least 4 systems is 

oo 3 

P(> 4, 1/) = ^ P(AA, i/) = 1 - ^ P(AA, i/) = 0.013 . (6) 

A/'=4 Af=0 

Similarly, if we take the number of systems observed towards SDSS J1148+5251 as a repre- 
sentative sample and compute the expected number of systems towards the remaining eight 
sightlines to he u — A {AX^f^^^/ AXff^g) ~ 23, then the hkelihood of detecting at most two 
systems is 

2 

P(< 2, 1/) = ^ P(j\^, i/) = 3.6 X 10-^ (7) 

Af=0 

The probability of simultaneously detecting at least 4 systems towards SDSS J1148-I-5251 
while detecting at most two systems along the remaining eight sightlines, 

2 



dX 



5:p(a/-,§^ax-3) 



Ar=o 



5:p(a/-,§^ax1j 



iAr=o 



(8) 



has a maximum of P{dM/dX) = 0.0018 at dU/dX = 0.84, taking a < 3. Thus, assuming 
that the O I systems are randomly distributed, we can rule out with a high degree of confi- 
dence the possibility that the excess number of detected systems towards SDSS J1148-I-5251 
is the result of a statistical fluctuation. 
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4.1. Compcirison with Lower- Redshift Sightlines 

Our sample of nine zqso > 4.8 sightlines is sufficient to demonstrate that the apparent 
excess of O I systems towards SDSS 1148+5251 is highly unlikely to arise from a random 
spatial distribution of absorbers. However, the distribution may be significantly non-random 
due to large-scale clustering. For example, if the line-of-sight towards SDSS 1148-1-5251 
passes through an overdense region of the IGM then we might expect to intersect more than 
the typical number of absorbers. 

To address this possibility, we have searched for O I in the spectra of an additional 30 
QSOs with 2.6 < zqso < 4.6. The sample was taken entirely with HIRES for a variety of 
studies. Due to the lower QSO redshifts, the data are of significantly higher quality than 
the Zqso > 4.8 sample. All sighthnes are complete down to A^o i — 10^^"^ cm~^ (apart from 
rare wavelength coverage gaps) with the majority complete down to A^oi = lO^^ '^ cm~^. 

A detailed study of the properties of lower-redshift O I absorbers will be presented 
elsewhere. Here we note that 11 systems with A^oi > 10^^'^ cm~^ were identified along all 
30 sightlines in the redshift intervals where O I would occur redward of the Lya forest. This 
is comparable to but slightly larger than the incidence rate among the eight lowest-redshift 
sightlines in our 2;qso > 4.8 sample. The number density is even somewhat higher among 
the -Zqso < 4.6 sample since the absorption pathlength intervals are shorter (cf. Eq. 2). 
However, many of the zqso < 4.6 spectra were taken in order to study known damped Lya 
systems (DLAs; Ahi > 10^°'^ cm~^), which commonly exhibit O I (e.g., Prochaska et al. 
2001, 2003). We would therefore expect these sighthnes to contain more than the typical 
number of O I systems. The important point is that no sightline contained more than a 
single system. If the excess towards SDSS J1148-I-5251 is due to the chance alignment with 
an overdense region of the IGM, then such a coincidence must occur in fewer than ~ 2 — 3% 
of all sighthnes. 



4.2. Comparison with Damped Lya and Lyman Limit System Populations 

We can compare the observed number of O I systems in our 2;qso > 4.8 sample with 
the expected number of DLAs by extrapolating from lower-redshift DLA surveys. Storrie- 

Lombardi & Wolfe (2000) find a line number density nDLA(-2) = dMnLA/dz = 0.055 {l+zY'^^. 
This is consistent at high- 2; with the larger survey done by Prochaska et al. (2005) and with 
the number density evolution Rao et al. (2005) obtained by combining DLA samples over 
< z < 4.5. From this t2dla(-2) we would expect a total of ~ 1.5 DLAs total along all nine 
lines-of-sight in the redshift intervals where are sensitive to O I. This is consistent with the 
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two O I systems seen along the eight lower-redshift sighthnes. However, we would expect 
to see only ~ 0.2 DLAs towards SDSS J1148+5251, which is inconsistent with the four we 
observe. 

The disparity eases for SDSS J1148+5251 if we include all Lyman hmit Systems (LLSs; 
A^Hi > 10^^-^ cm~^). The line number density of LLSs at 2; > 2.4 found by Peroux et 
al. (2003), nj^Lsiz) = 0.07(1 + z)^-'^^, predicts ~ 4 LLSs towards SDSS J1148+5251 in the 
interval where we can observe O I. However, there should also be ~ 20 LLSs along the 
remaining eight lines-of-sight. If the excess O I absorbers are associated with the progenitors 
of lower-redshift LLSs then they must evolve strongly in O I at 2 > 6. 



5. Metal Abundances 

Absorption systems at z > 5 provide an opportunity to study metal enrichment when the 
age of the Universe was < 1 Gyr. The opacity of the Lja forest at these redshifts prevents us 
from measuring reliable H I column densities. However, we are still able to constrain relative 
abundances. In Table 3 we summarize the Voigt profile fitting results. Total column densities 
for each system include only components with confirmed O I. Relative abundances were 
calculated using the solar values of Grevesse & Sauval (1998) and assuming zero ionization 
corrections. These are plotted in Figure 11. The error- weighted mean values for all six 
systems are ([0/Si]) = -0.04 ± 0.06, ([G/0]) = -0.31 ± 0.09, and ([G/Si]) = -0.34 ± 0.07 
(2(7 errors). 

The uncorrected abundances are broadly consistent with the expected yields from Type 
II supernovae of low-metaUicity progenitors (Woosley & Weaver 1995; Ghieffi & Limongi 
2004). Low- and intermediate-mass stars should not contribute significantly to the enrich- 
ment of these systems since these stars will not have had time to evolve by 2; ~ 6 (although 
see Barth et al. 2003). We can use the 0/Si ratio to constrain the enrichment from zero- 
metallicity very massive stars (VMSs, ~ 140 — 270 Mq). VMSs exploding as pair insta- 
bility supernovae are expected to yield ~ 4 times as much silicon compared to Type 11 SNe 
for the same amount of oxygen (Schaerer 2002; Umeda & Nomoto 2002; Heger & Woosley 
2002; Woosley & Weaver 1995; Ghieffi & Limongi 2004). We can estimate the fraction /^^^ 
of silicon contributed by VMSs to these systems using the simple mixing scheme of Qian & 
Wasserburg (2005a,b), 




Adopting [0/Si]vMS = -0.61 from Heger & Woosley (2002), and [0/Si]sNii = 0.00 from 
Woosley & Weaver (1995), the mean uncorrected 0/Si gives /^^^ = 0.11 ± 0.17 {2a). 
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However, this value of 0/Si should be considered a lower limit. The ionization potential of 
O I is smaller than that of Si ll (AEq j = 13.6 eV, AEg^ „ = 16.3 eV). Therefore, applying an 
ionization correction would increase 0/Si. Taking ([0/Si]) > —0.10 gives /^'^^ < 0.27 {2a). 

The mean C/Si ratio is consistent with zero silicon contribution from VMSs, although 
C/Si may decrease if ionization corrections are apphed {AEqu — 24:A eV > AE'sin). Dust 

corrections would also decrease C/Si since silicon is more readily depleted onto grains (Savage 
& Sembach 1996). Our upper limit of ([C/Si]) < —0.27 {2a) therefore does not provide an 
additional constraint on Z^^^^- However, the limit set by 0/Si may still imply that VMSs 
provided a smaller fraction of the metals in these O I systems than in the general IGM. Qian 
& Wasserburg (2005a, b) have argued that VMSs at z > 15 must have contributed > 50% of 
the sihcon in the IGM in order to produce [Si/C]iGM ~ 0.7 at 2; 2 — 4 (Aguirre et al. 2004). 
However, the Si/C ratio in the IGM derived from high- ionization hues depends sensitively 
on the choice of ionizing backgrounds. 

The high level of absorption blueward of the Lya emission hue at 2; > 5 prevents us 
from reliably fitting Voigt profiles to H I. However, transmission peaks in the Lyct and LyP 
forests near the O I redshifts allow us to set upper limits on the H I column density for 
two of our systems. Fits to these systems with the strongest allowable damping wings are 
shown in Figures 7 and 8, and the corresponding limits are listed in Table 3. The limits take 
into account uncertainties in the continuum level, which was set to a power-law F,y oc A^^ 
with spectral index (3 = 0.5, normalized at a rest wavelength of 1280 A. The continuum 
uncertainty will include the error in the amplitude of the power law at the wavelength of 
the damped system. Taking A/3 = 0.65 to incorporate 95% of the scatter in QSO spectral 
indices (Richards et al. 2001), this error is < 10% for the system at Zgys = 5.3364, for which 
the limit on A^h i was set using Lya transmission, and < 20% for the system at ^^gys = 6.0097, 
for which the limit was set using Ly/3. Larger uncertainties in the continuum may result from 
errors in the relative flux calibration and the departure from a power law. We conservatively 
estimate both of these to be ~ 20%, which translates to an uncertainty in the column density 
a{log Nui) ~ 0.3. The upper limits on log Ah i listed in Table 3 are our "best fit" values plus 
this uncertainty. The corresponding limits on the absolute abundances are [0/H] > —2.4 
(^,ys = 5.3364) and [0/H] > -4.1 {z^ys = 6.0097). 

The comoving mass density for each ion is given by 

where mion is the mass of the ion, pcrit = 0.92 x 10~^^ h^Q g cm~^ is the cosmological closure 
density, and AX*°* is the total pathlength interval over which we are sensitive to O I. In 
Table 4 we hst Qion for O l. Si 11, and C ll summing over all 2;qso > 4.8 sightlines, over the 
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three sightlines with zqso > 6.2 only, and over the sighthne towards SDSS J1148+5251 only. 
The mass densities include all confirmed absorption components, including those without 
detected O I. However, we have made no attempt to correct for incompleteness. The 
densities of low-ionization metal species along our zq^o > 6.2 sightlines are comparable to 
the densities of the corresponding high-ionization species at lower-redshift. Songaila (2005) 
measured il{C iv) ~ 3 x 10"*^ and fi(Si iv) ~ 1 x 10"^ over 1.5 < z < 4.5, with a drop in 
il{C iv) by at least a factor of two at z > 5. A decrease in high-ionization metals and a 
correlated rise in their low-ionization counterparts with redshift might signal a change in the 
ionization state of the z > 5 IGM. However, coverage of C IV does not currently extend to 
2; ~ 6, and our O I sample is probably too small to make a meaningful comparison. 



6. Discussion 

Although we observe an overabundance of O I at z > 6 along one of our highest- 
redshift sighthnes, it is not clear that these four systems constitute a "forest" which would 
signal a predominantly neutral IGM. In the reionization phase prior to the overlap of ionized 
bubbles. Oh (2002) predicts tens to hundreds of O I lines with A^oi > 10^^ cm~^, although 
the number may be lower if the enriched regions have very low metallicity {Z < 10'"^'^ Zq), 
or if the volume filling factor of metals is small {fz ^ 0.01). The number of O I systems 
wc detect is more consistent with his predictions for the post-overlap phase, although the 
expected number of lines depends on the ionizing efficiency of the metal-producing galaxies. 
The Oh (2002) model also predicts lines with A^^o i > 10^^ cm~^, which we do not observe. 

The significance of the excess O I towards SDSS Jl 148+5251 is complicated by the pres- 
ence of Lya and Ljf3 transmission along that sightlinc. The O I systems occur interspersed 
in redshift with the Ly/3 peaks, with the velocity separation from the nearest peak ranging 
from Av ~ 500 to 3000 km s~^ (AZ ~ 5 to 30 comoving Mpc). The O I absorbers also do 
not appear to be clustered. Along a total redshift interval A^tot = 0.49 (~ 200 comoving 
Mpc) where we are sensitive to O I, the detected systems are separated by Az — 0.06 to 0.12 
(AZ = 24 to 50 comoving Mpc). If any significantly neutral regions are present along this 
sightline, then they occur within or near regions that are highly ionized. In contrast, we find 
no O I towards SDSS J1030+0524. despite the fact that its spectrum shows complete ab- 
sorption in Lya and Ljj3 over Az ~ 0.2. The SDSS J1030+0524 sightline may pass through 
a region of the IGM that is still largely neutral and not yet enriched. Alternatively, there 
may be no O i because the sightline is highly ionized. In that case, a large-scale density 
enhancement could still produce the high Lya and Ly/5 optical depth (Wyithe & Loeb 2005). 

The overabundance of O I together with the transmission in Lya and Ly/3 indicates that 
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the SDSS J1148+5251 sightline has experienced both significant enrichment and ionization. 
One explanation is that we are looking along a filament. The total comoving distance 
between O I systems towards SDSS J1148+5251 is ~ 101 Mpc, which is smaller than the 
largest known filamentary structures at low redshift (e.g., Gott ct al. 2005). If the line- 
of-sight towards SDSS 1148+5251 runs along a filament at 2; ~ 6 then we might expect 
to intersect more than the average number of absorbers. Galaxies along such a structure 
might produce enough ionizing radiation to create the transmission gaps observed in the 
Lya and Ly/3 forests while still allowing pockets of neutral material to persist if the ionizing 
background is sufficiently low. If we assume that every 2; ~ 6 galaxy is surrounded by an 
enriched halo, then we can calculate the typical halo radius 

1/2 



-Rhalo 



7r(l + ^)2($g,i)5galADj 

35properkpc^^^"^^-^-V (^.0 ^ 



(11) 



4/ V 7.2 y VO-OlMpc-3 

\w) 1,200 Mpc J ■ 

Here, AD is the total comoving distance over which we are sensitive to O i. We have used 
the mean comoving number of galaxies ($gai) at 2; ~ 6 measured down to the limit of the 
Hubble Ultra Deep Field (UDF; Yan & Windhorst 2004; Bouwens et al. 2005), and assumed 
an overdensity of galaxies Sg^d — $gai/ (^gai) — 10 along this sightline. The absorption 
could be due to filled halos, shells, or other structures surrounding the galaxies. Note that 
by inverting Eq. 11 we can rule out the possibility that these O I systems arise from the 
central star-forming regions of galaxies similar those observed in the UDF. The small size 
-Rhaif-iight ~ 0.8 kpc of galaxics at 2; ~ 6 (Bouwens et al. 2005) would require either a local 
galaxy overdensity 5gai ~ 10'^ or a large number of faint galaxies to produce the observed 
absorption systems. 

At present wc can place relatively few physical constraints on these absorbers. Assuming 
charge- exchange equilibrium (Eq. 1), we can estimate the H I column densities as a function 
of metallicity, 

log iVH : = 16.9 + log (y^i^^^) - [0/H] . (13) 

For solar metallicities, even the weakest O I system in our sample must be nearly opti- 
cally thick. Metal-poor absorbers with [0/H] < —2 would require column densities A^hi ^ 
10^^ cm^^, which is consistent with the indication from the relative metal abundances that 
these absorbers are largely neutral. 



Finally, the excess O I absorbers are distinguished by their small velocity dispersions. 
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Three out of the four systems towards SDSS J1148+5251 have single components with 
Doppler widths 6 ~ 5 — 8 km s~^. We may be seeing only the low- ionization component of 
a kinematically more complex structure such as a galactic outflow (e.g., Simcoe et al. 2005). 
However, if there are additional high-ionization components that would appear in C IV and 
Si IV then we might expect to see components in C II and possibly Si II without O I. This 
is the case in the Zgys = 5.3364 system towards SDSS J0231-0728. The Zsys = 6.1293 and 
Zsys — 6.2550 systems towards SDSS J1148-I-5251 show potential C ll lines that do not appear 
in O I, but these are relatively weak, and in the Zgys — 6.1293 system there is no apparent 
strong Si IV. Among the 11 O I absorbers in the z < 4.6 sample, only three have a single 
O I component with a velocity width comparable to the three narrowest systems towards 
SDSS J1148-I-5251. The fact that the excess O I systems appear kinematically quiescent 
further suggests that they are a distinct class of absorbers. 



7. Summary 

We have conducted a search for O I in the high-resolution spectra of 9 QSOs at 4.8 < 
-^Qso < 6.4. In total, we detect six systems in the redshift intervals where O I A1302 falls 
redward of the Lya forest. Four of these lie towards SDSS Jl 148+5251 (zqso = 6.42). This 
imbalance is not easily explained by small number statics or by varying sensitivity among 
the data. A search at lower redshift revealed no more than than one O I system per sightline 
in a high-quality sample of 30 QSOs with 2.6 < zqso < 4.6. The number of excess systems is 
significantly larger than the expected number of damped Lya systems at 2; ~ 6, but similar 
to the predicted number of Lyman limit systems. 

The six systems have mean relative abundances ([0/Si]) = —0.04 ± 0.06, ([C/0]) = 
—0.31 ± 0.09, and ([C/Si]) = —0.34 ± 0.07 {2a errors), assuming no ionization correc- 
tions. These abundances are consistent with enrichment dominated by Type 11 supernovae. 
The lower limit on the oxygen to sihcon ratio ([0/Si]) > —0.10 is insensitive to ionization 
corrections and hmits the contribution of silicon from VMSs to < 30% in these systems. 
The upper limits on the ratio of carbon to oxygen [C/0] < —0.13, and carbon to silicon 
[C/Si] < —0.20 are also insensitive to ionization. Integrating over only the confirmed ab- 
sorption lines, the mean ionic comoving mass densities along the three zqso > 6.2 sightlines 
are Oqi ~ 7.0 x 10-^ Qsm ~ 9-6 x 10"^ and Qcu ~ 1-5 x 10"^ 

The excess O l systems towards SDSS Jl 148+5251 occur in redshift among Lya and LyP 
transmission features, which indicates that the IGM along this sightline has experienced both 
enrichment and significant ionization. The sightline may pass along a filament, where galax- 
ies are producing enough ionizing photons create the transmission gaps while still allowing a 
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few neutral absorbers to persist. In contrast, no O I was observed towards SDSS J1030+0524, 
whose spectrum shows complete Gunn- Peterson absorption. The SDSS J1030+0524 sightline 
may pass through an underdense region containing very few sources that would either enrich 
or ionize the IGM. Alternatively, the sightline may be highly ionized, with a large-scale den- 
sity enhancement producing strong line-blending in the Lja and Ljf3 forests. Future deep 
imaging of these fields and near-infrared spectroscopy covering high-ionization metal lines 
will allow us to distinguish between these scenarios and clarify the nature of the excess O I 
absorbers. 
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Table 1. Summary of Observations 



QSO 


ZQSO 


Dates 


Atexp (hrs)'' 


SDSS J2225-0014 


4.87 


2004 ,Jun 


6.7'' 


SDSS ,11204-0021 


5.09 


2005 ,Iaii - Feb 


6.7 


SDSS J0915+4244 


5.20 


2005 Jan - Feb 


10.8 


SDSS J0231-0728 


5.42 


2005 Jan - Feb 


10.0 


SDSS J0836+0054 


5.80 


2003 Feb - 2005 Jan 


21.7= 


SDSS J0002+2550 


5.82 


2005 Jan - Jun 


4.2 


SDSS J1623+3112 


6.22 


2005 Jun 


12.5 


SDSS J1030+0524 


6.30 


2005 Feb 


12.0 


SDSS J1148+5251 


6.42 


2003 Feb - 2005 Feb 


22.0"* 



''All exposures were taken with the upgraded HIRES detector unless 
otherwise noted. 

''Taken with the old HIRES detector 

"=9.2 hrs taken with the old HIRES detector 

"•4.2 hrs taken with the old HIRES detector 
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Table 2. Measured Properties of O I Systems 



Ion 


z 


b (km s-i) 


logiV (cm-2) 






SDSS .l(li>.ll — 07 




() 1 


Si II 


5.333982 ± 0.000008 


5.38 lb 0.58 


12.786 ±0.055 


P TTa 




u.Oo 




Si 11 


5.335928 ± 0.000020 


6.91 ± 1.33 


12.514 ± 0.079 


Si II 


5.336337 ± 0.000019 


7.62 ± 1.60 


12.978 ± 0.087 


Si ii'' 


5.336721 


12.96 


12.711 ±0.135 


O I'' 


5.335928 


6.91 


13.982 ± 0.086 


O 1=^ 


5.336337 


7.62 


14.138 ±0.129 


O I 


5.336721 ± 0.000072 


12.96 ±3.23 


13.972 ±0.143 


C 11^ 


5.335928 


6.91 


13.298 ± 0.077 


C ii« 


5.336337 


7.62 


13.454 ±0.119 


C ii'' 


5.336721 


12.96 


13.314 ±0.138 


Si II 


5.337841 ± 0.000014 


10.89 ±0.98 


12.678 ± 0.040 


Si II 


5.338474 ± 0.000043 


11.71 ±3.39 


12.233 ± 0.098 


Si II 


5.339301 ± 0.000016 


10.44 ± 1.13 


12.316 ± 0.063 


C 11^ 


5.337841 


10.89 


13.512 ± 0.051 


C 11^ 


5.338474 


11.71 


12.880 ±0.161 


C II* 


5.339301 


10.45 


13.074 ± 0.097 


Si W 


5.337841 


10.89 


12.848 ±0.051 


Si iv^ 


5.338474 


11.71 


12.142 ±0.192 


Si IV* 


5.339301 


10.44 


12.950 ± 0.049 


SDSS J1623+3112: Zsys = 5.8408 


O 1= 


5.836481 ±0.000012 


5.01 ±0.59 


14.376 ± 0.075 


O i'^ 


5.836806 ± 0.000017 


4.44 ± 1.57 


14.216 ±0.111 


Si li'^ 


5.836423 ± 0.000029 


3.63 ± 2.66 


12.746 ±0.182 


Si ii'^ 


5.836799 ± 0.000060 


10.82 ±4.29 


12.985 ±0.124 


Si 11^ 


5.839938 


6.37 


13.279 ±0.182 


Si 11^ 


5.840138 


4.13 


13.063 ±0.340 


Si ii'' 


5.840456 


6.92 


12.931 ±0.196 


Si lib 


5.840798 


5.32 


13.488 ± 0.060 


Si ll'' 


5.841113 


7.33 


13.271 ± 0.078 


Si 


5.841456 


1.77 


11.655 ± 1.018 


Si ii'* 


5.841780 


6.23 


12.854 ±0.070 


O 


5.839938 


6.37 


14.214 ±0.212 


O 


5.840138 


4.13 


14.152 ±0.242 


O I 


5.840456 ± 0.000032 


6.92 ±4.21 


13.734 ± 0.206 


O I 


5.840798 ± 0.000014 


5.32 ± 0.94 


14.384 ± 0.084 


O I 


5.841113 ± 0.000025 


7.33 ± 1.48 


14.316 ± 0.073 


O I 


5.841456 ± 0.000014 


1.77 ± 1.72 


13.518 ±0.150 


O i<i 


5.841780 


6.23 


13.793 ± 0.027 


C II 


5.839938 ± 0.000055 


6.37 ±1.56 


13.811 ±0.231 


C II 


5.840138 ± 0.000039 


4.13 ±2.50 


13.955 ± 0.249 
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Table 2 — Continued 



Ion 


z 


6 (km s-i) 


log AT (cm-2) 


C n'' 


5.840456 


6.92 


13.431 ±0.193 


C ii'' 


5.840798 


5.32 


13.771 ±0.093 


C ii'' 


5.841113 


7.33 


13.316 ±0.101 


C 11 


5.841780 ± 0.000007 


6.23 ±0.51 


13.736 ± 0.055 


SDSS J1U8+5251: Zsys = 6.0097 


Si ii'' 


6.009119 


7.02 


12.126 ±0.592 


Si li'' 


6.009290 


3.36 


12.770 ±0.137 


Si ll'' 


6.009717 


8.44 


13.323 ± 0.037 


Si lib 


6.010070 


4.94 


12.331 ±0.261 


O I 


6.009119 ± 0.000112 


7.02 ± 3.18 


13.593 ±0.371 


O I 


6.009290 ± 0.000012 


3.36 ±1.28 


14.071 ±0.113 


O I 


6.009717 ± 0.000010 


8.44 ± 0.96 


14.487 ± 0.030 


O I 


6.010070 ± 0.000034 


4.94 ± 1.60 


13.516 ±0.178 


C ii'' 


6.009119 


7.02 


13.180 ±0.328 


C ii'' 


6.009290 


3.36 


13.281 ±0.286 


C ii'' 


6.009717 


8.44 


14.013 ± 0.057 


C lib 


6.010070 


4.94 


13.044 ±0.182 


SDSS J1148+5251: Zsys = 6.1293 


C II': 


6.127673 ±0.000039 


6.11 ±3.03 


12.916 ±0.118 


C ii^: 


6.128589 ±0.000037 


8.08 ± 2.47 


13.101 ±0.092 


Si l& 


6.129267 


6.51 


13.273 ± 0.050 


O I 


6.129267 ±0.000006 


6.51 ±0.49 


14.606 ± 0.096 


C lib 


6.129267 


6.51 


13.744 ± 0.070 


C 11"= 


6.129885 ±0.000017 


2.21 ± 1.80 


13.161 ±0.186 




.s7),S'S' .iiii,f<+r, 




l<)(iS 


Si II 


6.196787 ±0.000014 


7.82 ± 0.84 


12.277 ± 0.044 


O i^ 


6.196787 


7.82 


13.727 ±0.061 


C 11^ 


6.196787 


7.82 


13.027 ± 0.072 


SDSS J1U8+5251: Zsys = 6.2555 


C II' 


0.255231 = O.OOOl.'W 


12.1::! ± 5.51 


1:1159 ± 0.229 


Si II 


6.255496 ± 0.000004 


4.30 ± 0.27 


12.900 ± 0.040 


O I 


6.255465 ± 0.000015 


5.08 ± 1.04 


14.075 ± 0.087 


C II* 


6.255496 


4.30 


13.596 ±0.114 



''The z- and 6-values for this component were tied to those for 
Si II. 

^The z- and 6-values for this component were tied to those for 
O I. 

Unconfirmed component 
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■^The z- and &- values for this component were tied to those for 
C II. 
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Table 3. Abundance Measurements 



Sightline 


^sys 


logATtota 


logJVtfi,- 


logJV*°*i- 


[O/Si]^ 


[C/O]^ 


[C/Si]^ 


log ATh I 


SDSS Jl 148+5251 


6.2555 


14.08 ± 0.09 


12.90 ±0.04 


13.60 ±0.11 


-0.11 ±0.08 


-0.17 ±0.15 


-0.27 ±0.12 




SDSS J1148+5251 


6.1968 


13.7.'? ±0.06 


12.28 ± 0.04 


13.03 ± 0.07 


0.17 ±0.07 


-0.39 ±0.10 


-0.22 ± 0.08 




SDSS J1148+5251 


6.1293 


14.61 ±0.10 


13.27 ±0.05 


13.74 ±0.07 


0.05 ± 0.10 


-0.55 ±0.13 


-0.50 ±0.09 




SDSS J1148+5251 


6.0097 


14.70 ±0.05 


13.48 ± 0.05 


14.17 ±0.07 


-0.07 ±0.06 


-0.22 ±0.09 


-0.28 ±0.08 


< 21.9= 


SDSS J1623+3112 


5.8408 


14.96 ± 0.06 


13.98 ± 0.06 


14.50 ± 0.09 


-0.31 ±0.09 


-0.15 ±0.13 


-0.45 ±0.11 




SDSS J0231-0728 


5.3364 


14.52 ±0.07 


13.25 ±0.06 


13.84 ± 0.07 


-0.02 ±0.07 


-0.37 ±0.12 


-0.39 ± 0.09 


< 20.1= 



^Total ionic column densities summed over all components with confirmed O I absorption 

''Relative abundances using the solar values of Grevesse & Sauval (1998) and assuming no ionization corrections 
"Upper limit including uncertainty in the continuum level 
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Table 4. Comoving Mass Densities' 



Sightlines 


n(o i) 


Q{Si II) 


Q{C II) 




(xlO-8) 


(xlO-9) 


(xlO-8) 


All 


3.2 ±0.2 


4.5 ±0.4 


0.8 ±0.1 


ZQSO > 6.2 


7.0 ±0.6 


9.6 ±0.9 


1.5 ±0.2 


SDSS J1148+5251 


10.2 ±1.0 


9.8 ±0.7 


1.8 ±0.2 



*For each ion, the value of Clion computed from Eq. 10 
includes all confirmed components listed in Table 2 along the 
indicated lines of sight, including those components without 
corresponding O I. These values arc lower limits in the sense 
that we have made no attempt to correct for incompleteness. 
Errors are la and reflect only the uncertainty in the summed 
column densities. 
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Fig. 1. — Absorption lines for the 2;sys = 5.3364 O I system towards SDSS J0231— 0728. 
Regions of the spectrum covering individual transitions have been shifted onto a common 
velocity scale. The histogram shows the normalized flux. The solid line shows the best Voigt 
profile fit. The dotted line shows the la flux error. 
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Fig. 2. — Absorption lines for the Zsys = 5.8408 O I system towards SDSS 1623+3112. See 
Figure 1 for details. The sharp features around Si ll A1260 are noise in the Lya forest. The 
components at A^; ~ —180 km s ^ are unconfirmed due to the strong Zabs — 

2.1979 Mg I 

absorption at the wavelength where we would expect to see C ll. 
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Fig. 3.— Absorption lines for the z^ys = 6.0097 O I system towards SDSS 1148+5251. See 
Figure 1 for details. Si II A1260 falls in the Lya forest. The strong line at A^; = 57 km s^^ 
from Si II A1304 is Si ll A1260 at z = 6.2555. The absorption features at Av ~ —40 km s~^ 
and Av pa 90 — 140 km s~^ from C ll are probably telluric absorption that has not been fully 
removed. 
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Fig. 4.— Absorption lines for the Zsys = 6.1293 O I system towards SDSS 1148+5251. See 
Figure 1 for details. The features around Si II A1260 are Lja absorption in the quasar 
proximity region. The C ll components at Av = —67, —29 and 26 km are unconfirmed. 
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Fig. 5. — Absorption lines for the Zgys = 6.1968 O I system towards SDSS 1148+5251. See 
Figure 1 for details. The features near Si II A 1260 are unrelated metal lines. The sharp 
feature near Si ll A1304 is the residual from a strong skyline. 



-31 - 



1.0 
0.5 
0.0 

1.0 
0.5 
0.0 

1.0 

0.5 

0.0 

1.0 
0.5 
0.0 

1.0 
0.5 
0.0 





1 






Si II 1260.42 ; 






Uj 










WW 


Si II 13 


D4.37 ; 


















V 






1 1302.17 : 
























V 




c 


Ir 

1 13 


^4 53 ; 


: n. H 1 1215.67 : 



100 

A V (km s~') 



100 



Fig. 6. — Absorption lines for the ^sys = 6.2555 O I system towards SDSS 1148+5251. See 
Figure 1 for details. The absorption at Av ^ —60 km s^^ from Si II A1260 is Si ll A1304 at 
z = 6.0097. The C ll fit includes an unconfirmed component on the blue edge of the profile. 
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Fig. 7. — Voigt profile fit to H I Lya and Ly/3 corresponding to the upper limit on A^h i for 
the ^sys = 5.3364 O I systems towards SDSS J0231— 0728. In order to set an upper hmit on 
A^H I we have taken b—1 km s~^, although the actual 6-parameter is likely to be significantly 
larger. In this system, A^h i is constrained by the transmitted flux in the wings of the Lya 
profile. 
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Fig. 8. — Voigt profile fit to H I Lya and Ly/3 corresponding to the upper limit on A^h i for 
the Zsys — 6.0097 O l systems towards SDSS J1148+5251. In this system, A^hi is constrained 
by the transmitted flux in the red wing of the Ly/3 profile. 
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Fig. 9. — Sensitivity curves for each of our sightlines. The solid lines give the absorption 
pathlength interval over which we are sensitive to single- component O I systems with A^^q i > 
Nq^^ and relative abundances similar to those of the O I systems in our sample. Circles mark 
the AX values corresponding to the A^oi of systems detected towards SDSS 1148+5251. 
The dashed line shows the mean AX as a function of Nq™ for all nine sightlines. The 
three lowest-redshift sightlines have smaller maximum AX values due in part to A-band 
atmospheric absorption. 
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Fig. 10. — Effective absorption pathlength interval over which we are sensitive to O I systems 
with > lO^'^'^ cm~^ as a function of a, where f{Noi) oc Nq'^ (see Eq. 3). The solid 
line shows AX^^ for SDSS J1148+5251. The dashed line shows the summed AX*^^ for the 
remaining eight sighthnes, divided by 5.7 to match SDSS J1148+5251 at a = 1.7. The dotted 
line shows the combined AX""^ for SDSS 1030+0524 and SDSS J1623+3112. For reasonable 
values of a, the ratio of AX''''^ for SDSS J1148+5251 and the remaining eight sightlines is 
essentially constant. We would therefore expect to see roughly six times as many O I systems 
along the other eight sightlines as towards SDSS J1148+5251 for a non-evolving population 
of absorbers. 
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Fig. 11. — Relative metal abundances in the O I systems assuming no ionization corrections. 
Individual measurements plotted as squares with la error bars are taken from Table 3. 
Dashed lines show the weighted mean values for all six systems. Dotted lines show the 2a 
uncertainties in the means. 



